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1 Introduction 

A new algebraic technique for solving Schrodinger and Klein-Gordon equations 
[1], and the related other works therein, has been introduced recently and used 
to search many interesting problems in different disciplines of physics. These 
works have clarified the power of the suggested model when compared the results 
obtained with those provided by other analytical methods in the literature. Nev- 
ertheless, this formalism involves a deficiency in its present form which requires, 
to express the excited state wave functions, the application of linear operators 
on the ground state wave function appeared automatically in the mathematical 
framework. This is indeed a cumbersome procedure though it provides explicit 
expressions for the state functions having non-zero angular momenta. 

To remove this drawback inherent in the formalism used in our previous works 
[1], we suggest here an alternative scheme, unifying the spirit of the two theoret- 
ical models [1, 2], to work out relativistic/non- relativistic quantum mechanical 
problems analytically in a unified framework. This is the main motivation be- 
hind the work presented in this article which in particular focuses on the solution 
of Dirac equation since recently considerable attention has been paid to exactly 
solvable Dirac equations. 

The arrangement of this article is as follows. In the next section, a brief 
introduction of the usual Dirac formalism and its treatment within the frame of 
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new scheme are presented. Third section involves apphcation results. Finally, 
the results are summarized in the concluding section. 



2 Theoretical Consideration 

To proceed, let us first remind the mathematical frame of Dirac formalism which 
is discussed briefly in the following section. Section 2.2 then illustrates the for- 
malism of the new model. 



2.1 Background on Dirac equation 

Dirac equation for scalar and vector potentials is given by [3] {h = c = 1) 

H^^{a.p + p{m + Vs) + Vv}^ (1) 

where p is the momentum operator, m is the rest mass of the particle, Vs and 
Vv are scalar and vector potentials respectively and a, (5 are Pauli matrices. 
To separate angular part of Eq. (1) from the radial part one follows 



iGti t 
— ^Ln. 
r T^jm 

L r r rjm J 



(2) 



where a represents the Pauli spin matrices while Gij and F^j are the radial com- 
ponents 

and (fi^j^ is the angular part of the wave function 

rjm 

Then, using the standard relations 



J =£-1/2 



(^■P-^V,m = --(^ + (3) 

and 

^ ^a.f gir) « i ,dg kg, , 



where 

U = ( + 1) = -0 + 1/2) i = f +1/2 1 ,^ ^ ^ ^^^^^^ 

\ +^ = +(j + l/2) J =£-1/2 J ^ 

Hence, the radial equations appear as 

dF{r) ^ k^^^^ = ^^_^_Vs- Vv)G{r) (5) 



dr r 
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+ -G(r) = {s^m + Vs- Vv)F{r), (6) 
ar r 

where e is the total relativistic energy of the system. From (5) and (6), omitting 
the derivatives of Vs and Vy, together with the ehmination of one radial compo- 
nent (-F(r)), one obtains a Schrodinger hke equation for the other component of 
the relativistic wave function 

{-^ + + i^s - V^) + {2mVs + 2eVy)^G = {e' - m')G. (7) 

2.2 Formalism 

Exact solutions of systems in physics has a great importance. To provide such 
solutions recently a new method has been carried out successfully in our earlier 
works [1], which unfortunately has a considerable algebraic difficulty in the cal- 
culation of excited state functions. To overcome this tedious procedure in the 
calculations, we propose here to use a new scheme involving orthogonal polyno- 
mials in order to express all bound state wave functions in an explicit form. 
For this purpose, we start from Eq.(7) which can be defined as 

G"{r) 

^ = VW-£, (8) 

where V{r) = + (\/| - V^) + {2mVs + leVy) and E = - m\ As is weU 

known, the solution of (8) generally takes the form 

G{r)^f{r)F[s{r)]. (9) 

The substitution of (9) into (8) yields the second-order differential equation 

and rearranging (10) for a more useful form, one gets 



Eq. (11) is in the form of the most familiar second-order differential equations 
to the hypergeometric type [2], 
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where a and a are at most second degree polynomials, and r is a first degree 
polynomial. The form of and is well defined for any special function 
F{s) [2]. From (12), it follows that 

s" t{s) f E-V a 

From the previous works in [1], the energy and potential terms in (13) can be 
decomposed in two pieces, which provides a clear understanding for the individual 
contributions of the F and / terms to the whole of the solutions, such that 
E - V = {Ef + Ep) - {Vf + Vf). Therefore, the second equality in (13) is 
transformed to a couple of equation 



ll^Vf-Ef, -^s''^Vf-Ef, (14) 



where / can be expressed in an explicit form due to the first part in (13) 

'1 t(s) 



/(r) = (s')-^/^exp 



2 J a(s) 



-ds 



(15) 



Since the corresponding a, a and r terms are well known for a given polynomial 
(F), the transformation function (s) in (14), and afterwards / in (15), are easily 
defined. So, from (9), the corresponding total wave function is readily obtained 
for the whole spectrum. 

The potential and total energy terms for Dirac equation in this case 



J-Vf- Ef, Vf = 2mVs + 2eVv + ^^^^J^ 



~s'''^Vf-Ef, Vf^V^-V^, (17) 



and 



Ef + EF = £^-m'^. {li 



To understand how efficiently this method works, some physically possible 
potentials are solved in the following section to obtain their eigenvalues and 
eigenf unctions, within the frame of the present formalism. 



3 Application 

As illustrative examples, here we deal with the two well known problems of the 
literature: Dirac oscillator and Dirac-Coulomb problem. The other solutions are 
shown in Table 1. 
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In order to get exact solutions for the present consideration, one needs to 
start with choosing a physicaUy plausible equal magnitudes for the vector and 
scalar potentials. Otherwise, the system considered becomes quasi-exactly solv- 
able which is out of the scope of the present article. 

Setting Vv = Vs — ar^ for the relativistic treatment of oscillator problem 
right hand sides of equations (16) and (17) gives 

y/ = 2a(m + £)r2 + M^±^ , Vp ^ , (19) 

Concentrating on the generalized Laguerre polynomials L'^{s) related to con- 
fluent hypergeometric functions, one sees that [2] 

a — s , T — a + 1 — s , a — na . (20) 



Keeping in mind that the right-hand sides of (14) provide a three dimensional har- 
monic oscillator potential, one obviously realizes that ^ = 2^ and consequently 
s — \ojr'^. Then, substituting (20) into (15) it is not hard to see that 



/ = Cr"+V2g-f^ ^ (21) 
where C — ;^(f )^- This makes possible to predict V) and as 

L. = V,-E,, Ef = (a+l)u,, 



where a = —{k + 1/2) = £ + | for the case j = £ + |. To find also Vp and Ep, 
one should consider Eq. (17). After some simple algebra we find 

Vf^O , EF^2nu . (23) 

Thus, in the non-relativistic limit, the full energy spectrum and wave functions 
for the system of interest are gives as 

E = Ef + Ef = {a + l + 2n)uj = (2n + £ + 3/2)w , 



* = /F = C.^e-iLr^\.). (24) 
Finally, the relativistic energy of Dirac oscillator reads 

£2 ^m2 + (2n + £ + 3/2)u;. (25) 
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The results obtained are in agreement with the work of Levai [2] which considers 
only the non-relativistic case, and also, for proper parameters {i + 1 — k, w — 
40:^^)5 these results agree well with the study of Alhaidari [4]. Moreover, the 
findings justify the excellent discussion in [5] on the confinement properties for 
Dirac equation with scalar and vector like potentials. 

Obviously, from the similarity between Eqs. (19) and (22) it is clear that 

" 8(m+£) ■ 

It is importantly noted that the choose of equal magnitudes for vector and 
scalar potentials leads to the non-relativistic limit of Dirac equation, removing 
the relativistic corrections. For a comprehensive understanding of this interesting 
point, the reader is referred to the individual works of Goniil and Kogak in [1] 
regarding the treatment of Klein-Gordon equation. We additionally remark that 
the present algebraic treatment has been performed only for spin-up case. Clearly, 
following similar procedure, one can easily repeat the same calculations for the 
other spinor wave function where now k — +i. 

The relativistic hydrogen atom is also an exactly solvable system within the 
frame of Dirac equation where the piece of potentials now should be defined as 

il^Vf-Ef, Vf^Vl-V^, (26) 



^Vf- Ef, Vf = 2mVs + 2eVv + M^±il . (27) 

g-2 J.Z 



For again the equal vector and scalar potentials; Vy = Vs = one gets 

In order to apply the present orthogonal polynomial technique, we choose the 
most suitable generalized Laguerre polynomial [F — e~*/^s~2~L^(s)] where 

2n + a + l l-a^ 1 
, = 1, , = 0. * = ^^ + ^-j, (29) 

2 

leading to s = ar, and to be in convenience with the literature, we set a = . 
Then, Eq. (15) reads 

f^t±l±2)i. (30) 
This justifies that — Vf — Ef — 0, while 

^^-"7 + ^^' ^^-"4(n + £ + l)2 ' 
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where k — —U+ 1) and a — 2i+l. Comparing Eqs. (28) and (31) one sees that 
b — 2(m+e) • Thus, for this system, the full energy spectrum and wave functions 
are 

E = Ef + E^ = E^ = -^^-^j—y^, (32) 



(33) 



where C — (ZLh^blH Finally, the relativistic energy for Dirac-Coulomb problem 
is 



4(n + £+l)2 ■ 



(34) 



The results are in agreement with the work of Levai [2] and, with a suitable 



aZ 
7+n+l 



1 -1 



overlap with those 



parameters such that 1 — 

obtained by Alhaidari in [4] . 

The formalism used here also provides explicit expressions for the relativistic 

spectra of three more potentials, which are illustrated in Table 1. We stress that, 
unlike the work in [4] , the strategy followed in this article for transforming Dirac 
equation into a Schrodinger like one does not enlarge the class of exactly solvable 
potentials as much as it might appear at first sight. Hence, the present brief work 
supports the individual criticism of the Castro and Vaidya-Rodrigues [6]. 

It is additionally stressed that one of the potentials listed in Table 1 do not 
require a restriction as in the examples discussed above in defining equal magni- 
tudes for vector and scalar potentials leading to the exact solvability. 



4 Concluding Remarks 

We have presented an idea of connecting the methods used in the analysis of ex- 
actly solvable potentials in non-relativistic quantum mechanics with the solution 
procedure of Dirac equation. The suggested formalism systematically recovers 
known results in a natural unified way and allows one to extend certain results 
known in particular cases. A straightforward generalization would be the ap- 
plication of the scheme to the other relativistic equations for integer spin cases. 
Beyond its intrinsic importance as a new solution for a fundamental equation in 
physics, we also expect that the present simple method would find a widespread 
application in the study of different quantum mechanical and nuclear scattering 
systems. Along this line the works are in progress. 
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Table 1: Relativistic energy and unnormalized eigenfunctions of the five poten- 
tials deduced within the present Dirac formalism discussed in section 3. In the 
treatment of Rosen-Morse and Eckart potentials, the notation carried out in [7] 
is used. 
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